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■ Abstract 

' Sub-Riemannian Geometry is proved to play an important role in many applications, e.g., Math- 

ematical Physics and Control Theory. Sub-Riemannian Geometry enjoys major differences from 

■ the Riemannian being a generalization of the latter at the same time, e.g., geodesies are not 
unique and may be singular, the Hausdorff dimension is larger than the manifold topological 
dimension. There exists a large amount of literature developing sub- Riemannian Geometry. How- 
ever, very few is known about its extension to pseudo-Riemannian analogues. It is natural to 
begin such a study with some low-dimensional manifolds. Based on ideas from sub-Riemannian 

■ geometry we develop sub-Lorentzian geometry over the classical 3-D anti-de Sitter space. Two 
different distributions of the tangent bundle of anti-de Sitter space yield two different geome- 
tries: sub-Lorentzian and sub-Riemannian. We use Lagrangian and Hamiltonian formalisms for 
both sub-Lorentzian and sub-Riemannian geometries to find geodesies. 
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II a ete prouve que la Geometrie sub-Riemannienne joue un role important dans des nom- 
breuses applications, par ex., Physique Mathematique et Theorie de Controle. La Geometrie sub- 
Riemannienne a des differences considerables par rapport a celle Riemannienne, etant au meme 
temps une generalisation de celle-ci, par ex., les geodesiques ne sont pas uniques et peuvent etre 
singulieres, la dimension de Hausdorff est plus grande que la dimension topologique de variete. 
II y a une quantite importante de litterature qui developpe la Geometrie sub-Riemannienne. 
Cependant, on connait tres peu sur son extension naturelle aux analogues pseudo Riemanniens. 
C'est naturel de commencer une telle etude avec des varietes de basse dimension. En se basant 
sur les idees de la geometrie sub-Riemannienne, on developpe la geometrie sub-Lorentzienne sur 
l'espace anti-de Sitter classique. Deux distributions differentes du faisceau tangent de l'espace 
d'anti-de Sitter donnent deux geometries differentes: sub-Lorentzienne et sub-Riemannienne. 
On utilise egalement les formalismes de Lagrange et d'Hamilton pour les deux geometries, sub- 
Lorentzienne et sub-Riemanniene, pour trouver les geodesiques. 
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1. Introduction 



Many interesting studies of anticommutative algebras and sub-Riemannian structures may be seen 
in a general setup of Clifford algebras and spin groups. Among others we distinguish the following 
example. The unit 3-dimensional sphere S 3 being embedded into the Euclidean space R 4 possesses a 
clear manifold structure with the Riemannian metric. It is interesting to consider the sphere S 3 as an 
algebraic object S 3 = SO(4)/SO(3) where the group SO (4) preserves the global Euclidean metric of 
the ambient space R 4 and SO(3) preserves the Riemannian metric on S 3 . The quotient SO(4)/SO(3) 
can be realized as the group SU(2) acting on S 3 as on the space of complex vectors z%, z-i of unit norm 
l z i| 2 + l z 2| 2 = 1- It is isomorphic to the group of unit quaternions with the group operation given by the 
quaternion multiplication. It is natural to make the correspondence between S 3 as a smooth manifold 
and S 3 as a Lie group acting on this manifold. The corresponding Lie algebra is given by left-invariant 
vector fields with non- vanishing commutators. This leads to construction of a sub-Riemannian structure 
on S 3 , see [4] (more about sub-Riemannian geometry see, for instance, |ll|19|20|21j ). The commutation 
relations for vector fields on the tangent bundle of S 3 come from the non-commutative multiplication for 
quaternions. Unit quaternions, acting by conjugation on vectors from R 3 (and R 4 ), define rotation in R 3 
(and R 4 ), thus preserving the positive-definite metric in R 4 . At the same time, the Clifford algebra over 
the vector space R 3 with the standard Euclidean metric gives rise to the spin group Spin(3) = SU(2) 
that acts on the group of unit spinors in the same fashion leaving some positive-definite quadratic form 
invariant. Two models are equivalent but the latter admits various generalizations. We are primary 
aimed at switching the Euclidean world to the Lorentzian one and sub-Riemannian geometry to sub- 
Lorentzian following a simple example similar to the above of a low-dimensional space that leads us 
to sub-Lorentzian geometry over the pseudohyperbolic space H 1 ' 2 in R 2,2 . In General Relativity the 
simply connected covering manifold of H 1,2 is called the universal anti-de Sitter space |15|16|22] , 

We start with some more rigorous explanations. A real Clifford algebra is associated with a vector 
space V equipped with a quadratic form Q(-, •). The multiplication (let us denote it by (g>) in the Clifford 
algebra satisfies the relation 

v (g) v = —Q(v, v)l, 

for v E V, where 1 is the unit element of the algebra. We restrict ourselves to V = R 3 with two different 
quadratic forms: 

r i 



Qs(v,v) = £v 



£ = 



and 



Q{v, v) = Iv ■ v, I = 



1 
1 

-10 
1 
1 



The first case represents the standard inner product in the Euclidean space R 3 . The second case corre- 
sponds to the Lorentzian metric in R 3 given by the diagonal metric tensor with the signature (—,+,+). 
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The corresponding Clifford algebras we denote by C1(0, 3) = Cl(3) and Cl(l, 2). The basis of the Clifford 
algebra Cl(3) consists of the elements 

{1, ii, i 2 , is, i\ <X> i 2 , h <$ H, i 2 <8> h, h®i 2 ® h}, with i\ ® ii = i 2 ® «2 = «3 ® h = — 1- 

The algebra Cl(3) can be associated with the space ixl, where H is the quaternion algebra. The basis 
of the Clifford algebra Cl(l, 2) is formed by 

{l,e,ii,»2, e® ii, e ® igjii ® *2>e® *i ® £2}) with e<8e = l, ii ® i\ = i 2 ® i 2 = — 1. 

In this case the algebra is represented by 2 x 2 complex matrices. 

Spin groups are generated by quadratic elements of Clifford algebras. We obtain the spin group 
Spin(3) in the case of the Clifford algebra Cl(3), and the group Spin(l,2) in the case of the Clifford 
algebra Cl(l, 2). The group Spin(3) is represented by the group SU(2) of unitary 2x2 complex matrices 
with determinant 1. The elements of SU(2) can be written as 

a,beC, H 2 + |6| 2 = l. 

The group Spin(3) = SU(2) forms a double cover of the group of rotations SO(3). In this case the 
Euclidean metric in M 3 is preserved under the actions of the group SO (3). The group Spin(3) = SU(2) 
acts on spinors similarly to how SO(3) acts on vectors from R 3 . Indeed, given an element R £ SO(3) 
the rotation is performed by the matrix multiplication RvR -1 , where «el 3 . An element U £ SU(2) 
acts over spinors regarded as 2 component vectors z = (zi,z 2 ) with complex entries in the same way 
Uzll^ 1 . This operation defines a 'half-rotation' and preserves the positive-definite metric for spinors. 
Restricting ourselves to spinors of length 1, we get the manifold {(zi, z 2 ) £ C 2 : \zi | 2 + \ z 2 \ 2 = 1} which 
is the unit sphere S 3 . 

Now we turn to the Lorentzian metric and to the Clifford algebra Cl(l, 2). The spin group Spin + (1, 2) 
is represented by the group SU + (1, 1) whose elements are 

a,kc, H 2 -|&| 2 = i. 

The group Spin + (1, 2) = SU + (1, 1) forms a double cover of the group of Lorentzian rotations SO(l, 2) 
preserving the Lorentzian metric Q(v, v). Acting on spinors, the group Spin + (1, 2) = SU + (1, 1) preserves 
the pseudo-Riemannian metric for spinors. Unit spinors (zi,z 2 ), \zi\ 2 — \z 2 \ 2 = 1, are invariant under 
the actions of the corresponding group Spin + (1,2) = SU + (1,1). The manifold if 1 ' 2 = {(zi,z 2 ) £ 
C 2 : \zi\ 2 — \z 2 \ 2 — 1} is a 3-dimensional Lorentzian manifold known as a pseudohyperbolic space 
in Geometry and as the anti-de Sitter space AdS^ in General Relativity. In fact, anti-de Sitter space 
is the maximally symmetric, simply connected, Lorentzian manifold of constant negative curvature. It 
is one of three maximally symmetric cosmological constant solutions to Einstein's field equation: de 
Sitter space with a positive cosmological constant A, anti-de Sitter space with a negative cosmological 
constant —A, and the flat space. Both de Sitter dSs and anti-de Sitter AdS$ spaces may be treated 
as non-compact hypersurfaces in the corresponding pseudo-Euclidean spaces R 1 ' 3 and M 2,2 . Sometimes 
de Sitter space dS^ or the hypersphere is used as a direct analogue to the sphere 5* 3 given its positive 
curvature. However, AdS^ geometrically is a natural object for us to work with. We reveal the analogy 
between S 3 and AdS3 as follows. The group of rotations SO (4) in the usual Euclidean 4-dimensional 
space acts as translations on the Euclidean sphere S 3 leaving it invariant. As it has been mentioned 
at the beginning, the sphere S 3 can be thought of as the Lie group S 3 = SO(4)/SO(3) endowed with 
the group law given by the multiplication of matrices from SU(2) which is the multiplication law for 
unit quaternions. The Lie algebra is identified with the left- invariant vector fields from the tangent 
space at the unity. The tangent bundle admits the natural sub-Riemannian structure and 5* 3 can be 
considered as a sub-Riemannian manifold. This geometric object was studied in details in [4]. It appears 
throughout celestial mechanics in works of Feynman and Vernon who described it in the language of 



a b 
—b a 



a b 
b a 
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two-level systems, in Berry's phase in quantum mechanics or in the Kustaaheimo-Stifel transformation 
for regularizing binary collision. 

Instead of R 4 , we consider now the space 

R 2,2 = { (xi, X2, X3, X4) G K 4 with a pseudo-metric dx 2 = —dx\ — dx\ + dx\ + dx 2 }. 

The group SO(2,2) acting on R 2 ' 2 is a direct analog of the rotation group SO(4) acting on R 4 . We 
consider AdS$ as a manifold H > = SO(2, 2)/ SO(l,2) with the Lorentzian metric induced from R 2,2 . 
Sometimes in physics literature, AdS^ appears as a universal cover of H > . It is worth to mention 
that H > is a homogeneous non-compact manifold and the group SO(2,2) acts as an isometry on 
H 1 ' 2 . The difference between this construction and above mentioned sphere is that S 3 itself is a group, 
whereas H 1 ' 2 is not. However, SO(2, 2) can be factorized as SO(2, 2) = SU+(1, 1) x SU+(1, 1)', so H 1 ' 2 
becomes a group manifold for SU + (1, 1), and topologically they are the same. The group law is defined 
by the matrix multiplication of elements from SU + (1, 1). The reader can find more information about 
the group actions and relation to General Relativity, e. g. [12ll7j . Left-invariant vector fields on the 
tangent bundle are not commutative and this gives us an opportunity to consider an analogue of sub- 
Riemannian geometry, that is called sub- Lorentzian geometry on SU + (1, 1) (which by abuse of notation, 
we call the AdS group). The geometry of anti-de Sitter space was studied in numerous works, see, for 
example, [H5ll0ll3ll8j . 

Very few is known about extension of sub-Riemannian geometry to its pseudo-Riemannian analogues. 
The simplest example of a sub-Riemannian structure is provided by the 3-D Heisenberg group. Let us 
mention that recently Grochowski studied its sub-Lorenzian analogue [7)8] . Our approach deals with 
non-nilpotent groups over 5 13 and AdS^. 

The paper is organized in the following way. In Section 2 we give the precise form of left-invariant 
vector fields defining sub-Lorentzian and sub-Riemannian structures on anti-de Sitter group. In Sections 
3 and 4 the question of existence of smooth horizontal curves in the sub-Lorentzian manifold is studied. 
The Lagrangian and Hamiltonian formalisms are applied to find sub-Lorentzian geodesies in Sections 
5 and 6. Section 7 is devoted to the study of a sub-Riemannian geometry defined by the distribution 
generated by spacelike vector fields of anti-de Sitter space. In both sub-Lorentzian and sub-Riemannian 
cases we find geodesies explicitly. 



2. Left-invariant vector fields 



We consider the AdS group topologically as a 3-dimensional manifold H 1 ' 2 in R 2 



H 1 ' 2 = {(xi,x 2 ,x 3 ,X4,) € 



oi = -i}, 



and as a group SU + (1, 1) with the group law given by the multiplication of the matrices from SU + (1, 1) 

a b 



We write a = x% + 1x2, b = X3 + 1x4, where i is the complex unity. For each matrix 



b a 



G SU + (1,1) 



we associate its coordinates to the complex vector p = (a, b). Then the multiplication law between 
p = (a, b) and q = (c, d) written in coordinates is 



pq = (a, 6)(c, d) = (ac + bd, ad + be). 



(2.1) 



Then, AdS with the multiplication law (|2.ip is the Lie group with the unity (1,0), with the inverse to 
p = (a, b) element p^ 1 = (a, —b), and with the left translation L p (q) — pq. The Lie algebra is associated 
with the left-invariant vector fields at the identity of the group. To calculate the real left-invariant vector 
fields, we write the multiplication law (12. ip in real coordinates, setting c = y± + iy%, d = j/3 + %4- Then 



pq = (xi, x 2 , x 3 , Xi)(yi, y 2 ,y3,y±) 

= (2:12/1 - x 2 y2 + x 3 y 3 + X4V4, x 2 yi + x x y 2 + x^yz - x 3 y 4 , 

£32/1 + £42/2 + £12/3 - £22/4, ^42/1 - £32/2 + ^22/3 + a:i2/4)- 



(2.2) 
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The tangent map corresponding to the left translation L p (q) is 



{Lp) t — 



xi -x 2 x 3 

X2 X\ X 4 

x 3 x 4 x\ 



x 4 

-x 3 

-X2 



X 4 



X 3 X 2 X\ 



The leftjinvariant vector fields are the left translations of vectors at the unity by the tangent map 
(i p )*: X = (L P )*X(0). Letting X(0) be the vectors of the standard basis in R 2 ' 2 (that coincides with 
the Euclidean basis in R 4 ), we get the left- invariant vector fields 

X\ = x\d Xl + x 2 d X2 + x 3 d X3 + x 4 d X4 , 
X 2 = -x 2 d Xl + xid X2 + x 4 d X3 - x 3 d Xi , 
X 3 = x 3 d Xl + x 4 d X2 + xid X3 + x 2 d X4 , 
X 4 = x 4 d Xl - x 3 d X2 - x 2 d X3 + x 1 d Xi 



in the basis d Xl ,d X27 d X3 ,d X4 . Let us introduce the matrices 



U 



10 
10 
10 
1 



10 
-10 
-1 
1 



E x = 



E 2 — 



1 
0-10 
0-1 
10 



10 
1 
10 
10 

Then the left-invariant vector fields can be written in the form 

X 1= xU-\7 X7 X 2 =xJ-\7 x , X 3 = xE 1 -V x , X 4 = xE 2 -\7 X} 

where x = (x\,x 2 ,X3,x 4 ), \7 X = (d Xl7 d X2 ,d X37 d Xi ) and "•" is the dot-product in R 4 . The matrices 
possess the following properties: 

• Anti-commutative rule or the Clifford algebra condition: 



JE 1 +E 1 J = 0, E 2 E 1 +E 1 E 2 ={), JE 2 + E 2 J = 0. 



Non-commutative rule: 



,11,1, ,1 r l 1,1 r l 1 , 1 

[^J,^E 1 }^-(JE 1 -E 1 J) = -E 2 , ^2 2 = 2 [ 2 J ' 2^ = " 2^' 



• Transpose matrices: 

• Square of matrices: 

As a consequence we obtain 

• Product of matrices: 



J — — J, E 2 — E 2 , E 1 — Ei. 



J 2 = -U, El = U, El = U. 
J E\ = E 2l E 2 Ei = J, JE 2 = — Ei. 



(2.3) 

(2.4) 

(2.5) 
(2.6) 

(2.7) 
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The inner (■, •) product in K 2,2 is given by 



(x, y) = Tx ■ y, with X = 




-10 
1 
1 



(2.8) 



Given the inner product |2 



we have 

(x,xEi) = (x,xJ) — {x,xE 2 ) =0, 
(xJ,xEi) — (xE2,xEi) — (xJ,xE 2 ) = 0, 
J,xJ)=-l, (xE 2 ,xE 2 ) = (xE 1 ,xEi) = 1. 



(2.9) 

(2.10) 
(2.11) 



The vector field X\ is orthogonal to AdS. Indeed, if we write AdS as a hypersurface F(x\ , x 2 , x 3 , X4 ) = 



1 = 0, then 



dF(c(s)) 



= 2 



dx\ dx 2 dx 3 dx4 
xi— x 2 — 1- x 3 — h x 4 



= (Xi, 



dc(s) N 



= 



ds \ ds ds ds ds I ' ds ' 

for any smooth curve c(s) = (xi(s), ^(s), X3(s), x^s)) on AdS. From now on we denote the vector field 
X\ by N. Observe, that |iV| 2 = (N, N) = —1. Up to certain ambiguity we use the same notation | • | as 
the norm (whose square is not necessary positive) of a vector and as the absolute value (non-negative) 
of a real/complex number. Other vector fields are orthogonal to N with respect to the inner product 
(•,-)inR 2 > 2 : 

(N,X 2 ) = (N,X 3 ) = (N,X 4 ) = 0. 

We conclude that the vector fields X 2l X3, X4 are tangent to AdS. Moreover, they are mutually 
orthogonal with 



IXol 2 = (Xn,Xo) = -1, 



\X 3 \ 2 = |X 4 | 2 = 1. 



We denote the vector field X 2 by T providing time orientation (for the terminology see the end of the 
present section). The spacelike vector fields X 3 and X4 will be denoted by X and Y respectively. We 
conclude that T, X, Y is the basis of the tangent bundle of AdS. In Table [T] the commutative relations 
between T, X, and Y are presented. We see that if we fix two of the vector fields, then they generate, 



Table 1 

Commutators of left-invariant vector fields 





T 


X 


Y 


T 





2Y 


-2X 


X 


-2Y 





-2T 


Y 


2X 


2T 






together with their commutators, the tangent bundle of the manifold AdS. 

Definition 1 Let M be a smooth n- dimensional manifold, T> be a smooth k- dimensional, k < n, 
bracket generating distribution on TM, and (-,-)x> be a smooth Lorentzian metric on T>. Then the 
triple (M,T>, (•, -)p) is called the sub- Lorentzian manifold. 

We deal with two following cases in Sections [3j-[6] and Section [7] respectively: 
1. The horizontal distribution V is generated by the vector fields T and X:T> — span{T, X}. In this case 
T provides the time orientation and X gives the spatial direction on D. The direction Y = ^[T,X], 
orthogonal to the distribution T>, is the second spatial direction on the tangent bundle. The metric 
(v)x> is given by the restriction of (•, •) from M 2,2 . This case corresponds to the sub-Lorentzian 
manifold {AdS,V, (-, -) v ). 



G 



2. The horizontal distribution T> is generated by the vector fields X and Y:T> — span{X, Y}. In this 
case both of the directions are spatial on T>. The direction T = i[Y, X], orthogonal to the distribution 
T>. In this case, the triple (AdS, T>, (•, -)-p) is a sub-Riemannian manifold. 

The ambient metric with the signature (—,—,+,+) of M 2,2 restricted to the tangent bundle TAdS of 
AdS is the Lorentzian metric with the signature (— , +, +), and therefore, AdS is a Lorentzian manifold. 
The vector fields T, X, Y form an orthonormal basis of each tangent space T p AdS at p £ AdS. We 
introduce a time orientation on AdS. A vector v £ T p AdS is said to be timelike if (v,v) < 0, spacelike 
if (v,v) > or v = 0, and lightlike if (v,v) = and ti ^ 0. By previous consideration we have T as a 
timelike vector field and X, Y as spacelike vector fields at each p £ AdS*. A timelike vector v £ T p AdS 
is said to be future-directed if (w, T) < or past-directed if (v, T) > 0. A smooth curve 7 : [0, 1] — > AdS 1 
with 7(0) = p and 7(1) = g is called timelike (spacelike, lightlike) if the tangent vector j(t) is timelike 
(spacelike, lightlike) for any t £ [0, 1]. If £l Pt<] is the non-empty set of all timelike, future-directed smooth 
curves j(t) connecting the points p and q on AdS, then the distance between p and q is defined as 

1 

dist := sup J y/-(j{t),j(t))dt. 


A geodesic in any manifold M is a curve 7 : [0, 1] —> M whose vector field is parallel, or equivalently, 
geodesies are the curves of acceleration zero. A manifold M is called geodesically connected if, given 
two points p,q £ M, there is a geodesic curve j(t) connecting them. Anti-de Sitter space AdS is not 
geodesically connected, see [9114] . 

The concept of causality is important in the study of Lorentz manifolds. We say that p £ M chrono- 
logically (causally) precedes q £ M if there is a timelike (non-spacelike) future-directed (if non-zero) 
curve starting at p and ending at q. For each p £ M we define the chronological future of p as 

I + {p) = {q £ M : p chronologically precedes q}, 

and the causal future of p as 

J + (p) = {q £ M : p causally precedes q}. 

The conformal infinity due to Penrose is timelike. One can make analogous definitions replacing 'future' 
by 'past'. 

From the mathematical point of view the spacelike curves have the same right to be studied as timelike 
or lightlike curves. Nevertheless, the timelike curves and lightlike curves possess an additional physical 
meaning as the following example shows. 

Example 1. Interpreting the ^-coordinate of AdS as time measured in some inertial frame (xi = t), 
the timelike curves represent motions of particles such that 

(£)' + (£)'<•• 

It is assumed that units have been chosen so that 1 is the maximal allowed velocity for a matter particle 
(the speed of light). Therefore, timelike curves represents motions of matter particles. Timelike geodesies 
represent motions with constant speed. In addition, the length 

1 

r(l) = J y/-W),i(t))dt, 


of a timelike curve 7 : [0, 1] — > AdS is interpreted as the proper time measured by a particle between 
events 7(0) and 7(1). 

Lightlike curves represent motions at the speed of light and the lightlikc geodesies represent motions 
along the light rays. 
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3. Horizontal curves with respect to the distribution T> — span{T, X} 



Up to Section [7] we shall work with the horizontal distribution T> = span{T, X} and the Lorentzian 
metric on V, which is the restriction of the metric (•, •) from R 2 ' 2 . We say that an absolutely continuous 
curve c(s) : [0, 1] —* AdS is horizontal if the tangent vector c(s) satisfies the relation c(s) = a(s)T(c(s)) + 
P(s)X(c(s)). 

Lemma 1 A curve c(s) = (xi(s), X2(s), x$(s), Xi(s)) is horizontal with respect to the distribution T> = 
span{T, X}, if and only if, 

— X4X1 + X3X2 — X2X3 + X1X4 = or (xE2,c) = 0. (3-1) 

PROOF. The tangent vector to the curve c(s) = (xi(s), 0:2 (s), ^(s), Xi{a)) written in the left-invariant 
basis (T, X, Y) admits the form 

c(s) = ctT + (3X + jY. 

Then 

7 = (c, Y) = Xc ■ Y = -.X4.i1 + x 3 x 2 — X2X3 + x\Xi = (xE 2 , c). 

We conclude that 

7 = 0, 

if and only if, the condition (|3.1|) holds. □ 

In other words, a curve c(s) is horizontal, if and only if, its velocity vector c(s) is orthogonal to 
the missing direction Y. The left-invariant coordinates ct(s) and /3(s) of a horizontal curve c(s) = 
(xi(s),X2(s),x 3 (s),x i (s)) are 

a = (c, T) = X2X1 — x\X2 + X4X3 — X3X4 — (xJ, c), (3-2) 

(3 = (c, X) = —X3X1 — X4X2 + X1&3 + X2X4 — (xE\, c). (3-3) 

Let us write the definition of the horizontal distribution T> = span{T, X} using the contact form. We 
define the form u = —x^dxi + x$dx2 — £2^3 + x\dx^ = (XE2, dx). Then, 

oj{N) = 0, lu(T) = 0, uj(X) = 0, lu(Y) = 1, 

and kerw = spanjiV, T, Y}, The horizontal distribution can be defined as follows 

V = {V € TAdS : lu{V) = 0}, or V — ker oj n TAdS, 

where TAdS is the tangent bundle of AdS. 

The length 1(c) of a horizontal curve c(s) : [0, 1] — > AdS is defined by the following formula 

i(c) = f'msicwi^ds. 

Jo 

Using the orthonormality of the vector fields T and X, we deduce that 

1(c) = [ \~a 2 (s)+f3 2 (s)\ 1/2 ds. 
Jo 

We see that the restriction onto the horizontal distribution T> C TAdS of the non-degenerate metric 
(•, ■) defined on TAdS gives the Lorentzian metric which is non-degenerate. The definitions of timelike 
(spacelike, lightlike) horizontal vectors v £ V p are the same as for the vectors v € T p AdS. A horizontal 
curve c(s) is timelike (spacelike, lightlike) if its velocity vector c(s) is horizontal timelike (spacelike, 
lightlike) vector at each point of this curve. 

Lemma 2 Letj(s) — (yi(s),y2(s),ys(s),y4,(s)) be a horizontal timelike future- directed (or past- directed) 
curve and c(s) = L p ("f(s)) be its left translation by p = (pi,P2,P3,Pi), P € AdS. Then the curve c(s) is 
horizontal timelike and future- directed (or past-directed). 
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PROOF. Let us denote by (ci(s), e^s), C3(s), Ci(s)) the coordinates of the curve c(s). Then, by (|2.2|) 
we have 



Differentiating with respect to s, we calculate the horizontality condition (|3.ip for the curve c(s). Since 
— p\ — p\ + v\ + p\ — — 1, straightforward simplifications lead to the relation 

(c, K) = -C4C1 + c 3 c 2 - c 2 c 3 + cic 4 = (-p^ - pi +pj+pl)(-y4yi + 2/32/2 - 2/22/3 + 2/12/4) = 0, 

and the curve 7 is horizontal. 

Let us show that the curve c(s) is timelike and future-directed provided 7(5) is such. We calculate 

(c,T) = c 2 ci - cic 2 + C4C3 - C3C4 = (pi +p\ -p\ -pi) (2/22/1 - 2/12/2 + 2/42/3 - 2/32/4) = (i,T) 



(c,X) = -C3C1 - c 4 c 2 +C1C3 + c 2 c 4 = (p? +^2 -p§ (-2/32/1 - 2/42/2 + 2/12/3 + 2/22/4) = (7, A) 

from (|3.2[) . (|3.3p . and (|3 .4|) . Since the horizontal coordinates are not changed, we conclude that the 
property timelikeness and future-directness is preserved under the left translations. □ 

In view that the left-invariant coordinates of the velocity vector to a horizontal curve do not change 
under left translations, we conclude the following analogue of the preceding lemma. 

Lemma 3 Let^{s) = (yi(s), 2/2(s), 2/3(5), 2/4( s )) ^ e a horizontal spacelike (or lightlike) curve and c(s) = 
Lp("f(s)) be its left translation by p = (j>i,P2,P3,P4), P S AdS . Then the curve c(s) is horizontal spacelike 
(or lightlike). 

4. Existence of smooth horizontal curves on AdS 

The question of the connectivity by geodesies of two arbitrary points on a Lorentzian manifold is 
not trivial, because we have to distinguish timelike and spacelike curves. The problem becomes more 
difficult if we study connectivity for sub-Lorentzian geometry. In the classical Riemannian geometry all 
geodesies can be found as solutions to the Euler-Lagrange equations and they coincide with the solutions 
to the corresponding Hamiltonian system obtained by the Legendre transform. In the sub-Riemannian 
geometry, any solution to the Hamiltonian system is a horizontal curve and satisfies the Euler-Lagrange 
equations. However, a solution to the Euler-Lagrange equations is a solution to the Hamiltonian system 
only if it is horizontal. 

In the case of sub-Lorentzian geometry we have no information about such a correspondence. As it 
will be shown in Sections 6 and 7 the solutions to the Hamiltonian system are horizontal. It is a rather 
expectable fact given the corresponding analysis of sub-Riemannian structures, e. g., on nilpotent groups, 
see |2I3) . Since {T, X, Y = 1/2[T,X]} span the tangent space at each point of AdS the existence of 
horizontal curves is guaranteed by Chow's theorem [6]. So as the first step, in this section we study 
connectivity by smooth horizontal curves. The main results states that any two points can be connected 
by a smooth horizontal curve. A naturally arisen question is whether the found horizontal curve is 
timelike (spacelike, lightlike)? 

First, we introduce a parametrization of AdS and present the horizontality condition and the hori- 
zontal coordinates in terms of this parametrisation. 

The manifold AdS can be parametrized by 



c 2 (s) 
cs(s) 
c 4 (s) 



Pl2/l0) -P22/2(s) +P32/3(«) +P42/4(s), 
P22/l(s) +P12/20) +P42/3(s) -P3Vi[s), 
P3Vl(s) +P42/2(s) +Pm{s) -P2Va{s), 
P4Vl(s) -P3V2(s) +P22/3(s) +PiVa(s)- 



(3.4) 



and 
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X\ = cos a cosh 0, 
X2 = sin a cosh 0, 

(4.1) 

xs = cos&sinh0, 
X4, = sm&sinh0, 

with a, b G (— 7r, +7r], € (— oo, oo). Setting ip = a — b, ip = a + b, we formulate the following lemma. 

Lemma 4 Let c(s) = (f(s), il>(s), 6(s)) be a curve on AdS. The curve is horizontal, if and only if, 

V9cos-0sinh26» - 20 sin -0 = 0. (4.2) 
The horizontal coordinates a and (3 of the velocity vector are 

a = --(<pcosh20 + ^>) = -dcosh 2 9 - 6sinh 2 9, (4.3) 
(3 = -(<p sin ^sinh20 + 29cosip). (4.4) 

PROOF. Using the parametrisation (|4.ip of AdS, we calculate 

i\ = — d sin a cosh + 0cosasinh0, 

Xi = d cos a cosh + 0sinasinh0, 

±3 = — &sinfesinh0 + cos & cosh 0, 

X4 = &cos&sinh0 + sin 6 cosh 0. 

Substituting the expressions for Xk and ik, k= 1,2, 3, 4, in (|3.ip . (|3.2p . and (|3.3p . in terms of ^ and 
9, we get the necessary result. □ 

We also need the following obvious technical lemma formulated without proof. 
Lemma 5 Given qo, q\, Jel, there is a smooth function q : [0, 1] — > K ; smc/i t/iat 

9(0) =q , ?(l)=?i, / q(u)du = I. 

Jo 

Theorem 1 Let P cmrf Q be two arbitrary points in AdS. Then there is a smooth horizontal curve 
joining P and Q. 

PROOF. Let P = P((fo, V'o, #o) and Q — Q((fi, tpi, 0i) be coordinates of the points P and Q. In order 
to find a horizontal curve c(s) we must solve equation (|4.2p with the boundary conditions 

c(0) = P, or <^(0)=^o, ^(0)=^o, 0(0) = 9 , 

c(l) = Q, or tp{l)=<pi, V(l)=^i, 0(1) = ^i. 

Assume that sin?/; ^ we rewrite the equation (|4. 2[) as 

20 = (^coti/isinh20. (4.6) 

To simplify matters, let us introduce two new smooth functions p(s) and q(s) by 

20(s) = arcsinhp(s), t/j(s) — arccot q(s), 

and let the function ip(s) is set as y(s) = ifo + s(tpi — ipg). Then we will define the smooth functions p(s) 
and q(s) satisfying the horizontality condition (|4.6p for c = c(s). Let k = tpi — ipo. Then equation (14. 6p 
admits the form 

7f=f7 \ = k P( s M s )- 
i/l +p 2 (s) 
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Separation of variables leads to the equation 

dp 



kq(s) ds, 



p^l+p 2 

that after integrating gives 

— arctanh — - 

To define the constant C, we use the boundary conditions at s = 0. Observe that 

and 



fc( f q{T)dr + Cy 



v/l+p 2 (0) cosh20 o v/l+p 2 (l) cosh26»i 

Then 

1 1 
C — — — arctanh - 



A: cosh 20o 

Applying the boundary condition at s = 1 we find the value of f Q q(r) dr as 

/ q(t) dr = ( arctanh h arctanh 

J k\ cosh 20i cosh 26»o 

Since, moreover, q(0) = coti/>o, <z(l) = coti/'i, Lemma [5] implies the existence of a smooth function g(s) 
satisfying the above relation. 

The function p(s) can be defined by 



VTTpW) 



tanh 



f 1 

k / ofr) dr — arctanh — „ 

Jo cosh 20 o 



The curve c(s) = (tp(s), -0(s), 0(s)) = (<y5o + s(^i — ipo)> arccot q(s)), \ arcsinhp(s)) is the desired 
horizontal curve. □ 



Remark 1 Of course, the proof is given for a particular parametrisation by a linear function Lp. One 
may easily modify this proof for an arbitrary smooth function ip obtaining a wider class of smooth 
horizontal curves. 

Some of the points on AdS can be connected by a curve that maintain one of the coordinate constant. 



Theorem 2 If P = P(ipo,ip,0o) and Q = Q((fi, if), 0i) with 

r tanh 0i , , \ 

^ arCC0 U ln ^nW^°-^)J (4J) 

are two points that can be connected, then there is a smooth horizontal curve joining P and Q with the 
constant ip- coordinate given by (|4.7p . 



PROOF. Let c = c(cp, ip, 8) be a horizontal curve with the constant ^-coordinate. Then it satisfies the 
equation (|4.2p that in this case we write as 

, , d(20) 
cot ip dip — — 



sinh 20 
Integrating yields 

rd d(26) 



cot ip / dip 

100 



sinh 20 

cot ip((p(6) - <p(0 )) = In tanh - In tanh O . (4.8) 
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For 9 = 6\ we get formula (|4.7p for the value of ip. Solving (|4. 8|) with respect to <p(9) we get 

In ( tanh 9/ tanh On) 

with ip given by (|4.7p . Finally, the horizontal curve joining the points P and Q satisfies the equation 



/ In ( tanh 6/ tanh An) \ 



cot ^ 

□ 

Upon solving the problem of the connectivity of two arbitrary points by a horizontal curve we are 
interested in determining its character: timelikeness (spacelikeness or lightlikcness) . It is not an easy 
problem. We are able to present some particular examples showing its complexity. Let us start with the 
following remark. 

Remark 2 If P,Q € AdS are two points connectable only by a family of smooth timelike ( spacelike, 
lightlike) curves, then smooth horizontal curves (its existence is known by the preceding theorem) joining 
P and Q are timelike (spacelike, lightlike). 

Indeed, let flp,Q be a family of smooth timelike (lightlike) curves connecting P and Q. If S(s) S Qp,Q, 
then its velocity vector S(s) can be written in the left-invariant basis T, X, Y as 

5(s) = a(s)T(S( S ))+f3(s)X(S(s))+ 1 (s)Y(8(s)) 

with (8(s),8(s)) = —a 2 + (3 2 + j 2 < 0(= 0). If moreover, it is horizontal, then 7 = 0. Therefore, 
—a 2 + (3 2 < 0(= 0), and the horizontal curve connecting P and Q is timelike (lightlike). 

If the points P and Q are connectable only by a family of spacelike curves, then the inequality 
—a 2 +(3 2 > 7 2 holds for them. It implies — a 2 + (3 2 > for a horizontal curve. We conclude that in this 
case the horizontal curve is still spacelike. 

Making use of (|4.3|) and (|4.4p as well as parametrisation (|4.1j) we calculate the square of the velocity 
vector for a horizontal curve in terms of the variables tp, ip, 9 as 

-a 2 +/3 2 = -ip 2 - + A9 2 - 2<pip cosh 29. (4.9) 
We present some particular timclikc, spacelike, and lightlike solutions of (|4.2[) . 



Example 2. Let tp = 0. Then, tp = <p is constant. In order to satisfy (|4.2[) we have two options: 

2.1 9 = => 9 = 9q is constant. Then |c| 2 = — ip 2 < 0. We conclude that all non-constant 
horizontal curves c(s) = (tpo, tp( s )i &o) are timelike. The projections of these curves onto the {x\, X2)- 
and (X3, X4) -planes are circles. All lightlike horizontal curves are only constant ones. 

2.2 ip = nn, neZ. Then |c| 2 = 4(9 2 > 0. We conclude that all non-constant horizontal curves c(s) = 
(tpQ,nn,9(s)), n S Z are spacelike. The projections of these curves onto the (0:1,2:3)- and (2:2,2:4)- 
planes are hyperbolas. All lightlike horizontal curves are only constant ones. 



Example 3. Let tp ^ 0. We choose tp as a parameter. Then the square of the norm of the velocity 
vector is 

- a 2 + I3 2 = -1 - ^P 2 + 4<9 2 - 2<ip cosh 29, (4.10) 
where the derivatives are taken with respect to the parameter tp. The horizontality condition becomes 

26>sin^> = cos^smh26>. (4.11) 

As in the previous example we consider different cases. 
3.1 Suppose 9 = and assume that 9 = 9q 0. Then the horizontal curves are parametrized by 
c(s) = (tp, I + 7m, #0)7 n 6 Z. All these curves are timelike, since |c| 2 = — 1. There are no lightlike or 
spacelike horizontal curves. 
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3.2 If #o = 0, then any curve in the (tp, ^/>)-plane is horizontal and timelike since |c| 2 = —(1 + ip) 2 - 

3.3 Suppose that ip = and i/> e i/> ^ f , fc e Z. Then (|4~T0f and (|4~TT]) are simplified to 

-a 2 +p 2 = -1 + 49 2 , (4.12) 

8 = Ksmh26 with K= ^ll ( 4 13 ) 
Let = 9(ip) solves equation (|4. 1 3|) . Then the horizontal curve 

c(s) = (^ ,%)) (4.14) 

is timelike when \8\ < iarcsinh-^. If \9\ > (=)i arcsinh then the horizontal curve (14.14[) is 
spacelike (lightlike). 

Thus any two points P(ip Q ,tpQ,9o), Q( l Pii'4 > ii@o)j can be connected by a piecewise smooth timelike 
horizontal curve. This curve consists of straight segments with constant (^-coordinates or with coordinate 
if) = ? + im, n € Z. In the case 9q — 0, this horizontal curve can be constructed to be smooth. 

5. Sub-Lorentzian geodesies 

In Lorentzian geometry there are no curves of minimal length because two arbitrary points can be 
connected by a piecewise lightlike curve. However, there do exist timelike curves with maximal length 
which are timelike geodesies [14] . By this reason, we are looking for the longest curve among all horizontal 
timelike ones. It will be given by extremizing the action integral S = \ J Q ( — a 2 (s) + /3 2 (s)) ds under 
the non-holonomic constrain (xE2, c) = 0. The extremal curve will satisfy the Euler-Lagrange system 

d dL _ dL ^ 
ds dc dc 

with the Lagrangian 

L(c,c) = i(~a 2 +P 2 ) + \{s){xE 2 ,c). 

The function A(s) is the Lagrange multiplier function and the values of a and are given by (13. 2|) 
and (|3.3p . The Euler-Lagrange system (|5.ip can be written in the form 

— &X2 — (3x3 = 2{a±2 + Px^ — XX4) — Xx4, 

axi — I3xa = 2(— axi + (3x4 + ^±3) + Xx^, 
—61x4 + $xi = 2{a±4 — (3x\ — Xx,2) — Aa;2, 

ax^ + (3x2 = 2(— ais — (3x2 + Xi\) + XX4. 

for the extremal curve c(s) = (xi(s),X2(s),X3(s),Xi(s)). Multiplying these equations by X2, —Xi, —X4, 
X3, respectively and then, summing them up we obtain 

-a = 2(-a(c, N) - (3{c, Y) - A/3) = -2A/3 

because (c, Y) — (c,N) — 0. Now, multiplying the equations by X3, X4, x\, x 2 , respectively and then, 
summing them up we get 

-(3 = 2(a(c, Y) + (3{c, N) + Xa) = 2Xa 
in a similar way. The values of a and (3 are concluded to satisfy the system 

a(s) = 2X(3(s) 1 

W W (5.2) 

/3(s) = 2Aa(s). 

Case X(s) = 0. In the Riemannian geometry the Schwartz inequality allows us to define the angle t9 
between two vectors v and w as a unique number < i3 < it, such that 

v ■ w 

COS 17 = 7— 7| r . 

\v\\w\ 

There is an analogous result in Lorentzian geometry which is formulated as follows. 
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Proposition 1 J l^j Let v and w be timelike vectors. Then, 

1. \{v, w}\ > \v\\w\ where the equality is attained if and only if v and w are collinear. 

2. If (v, w) < 0, there is a unique number $ > 0, called the hyperbolic angle between v and w, such that 

(v, w) = — \v\\w\ coshi?. 

Theorem 3 The family of timelike future- directed horizontal curves contains horizontal timelike future- 
directed geodesies c(s) with the following properties 

1. The length \c\ is constant along the geodesic. 

2. The inner products (T, c) = a, (X, c) — (3, (Y, c) = are constant along the geodesic. 

3. The hyperbolic angle between the horizontal time vector field T and the velocity vector c is constant. 

PROOF. The system (jO} implies 

a(s) = $(s)=0. 

The existence of a geodesic follows from the general theory of ordinary differential equations, employ- 
ing, for example, the parametrisation given for a, 0, 7 in the preceding section. Since the horizontal 
coordinates a(s) and (3(s) are constant along the curve c we conclude that c is geodesic. We denote by 
a and f3 its respective horizontal coordinates. 

The length of the velocity vector c is | — a 2 + ft 2 ] 1 / 2 and it is constant along the geodesic. 

The second statement is obvious. Since c(s) is a future-directed geodesic, we have (T, c) < 0, and 

w , x (T,c) -a 

coshlZT, c) = — }—. , = — is constant. 

m|c| ^\-a 2 +f3 2 \ 

□ 

Case A(s) 7^ 0. We continue to study the extremals given by the solutions of the Euler-Lagrange 
equation (|5.ip . 

Lemma 6 Letc(s) be a timelike future- directed solution of the Euler-Lagrange system (|5.ip with A(s) ^ 

0. Then, 

1. The length \ — a 2 (s) + P 2 (s)\ 1 ^ 2 of the velocity vector c(s) is constant along the solution. 

2. The hyperbolic angle between the curve c(s) and the integral curve of the time vector field T is given 
by 

§ = Z(c,T) = -2A(s) + O , 

where A is the primitive of A. 

PROOF. Multiplying the first equation of (|5.2p by a, the second one by f3 and subtracting, we deduce 
that aa — /30 = 0. This implies that —a 2 + (3 2 = (c, c) is constant. The horizontal solution is timelike 
if the initial velocity vector is timelike. The first assertion is proved. 
Set r = y| ~ 0/2 + P 2 \- Using the hyperbolic functions we write 

a(s) — —r cosh 9(s), /3(s) = rsinh#(s). 

Substituting a and fi in (|5.2j) . we have 

9(s) = -2A(s). 

Denote A(s) = A(s) ds and write the solution of the latter equation as 9 = — 2A(s) + 9q. Thus, 

a(s) = -rcosh(-2A(s) + O ), /9(s) = r sinh(-2A(s) + 9 ). (5.3) 

In order to find the value of the constant 9q we put s — and get 9 = arctanh . 
Let c(s) be a horizontal timelike future-directed solution of (|5.ip . Then (c, T) < and 

a = (c,T) = -|c||T|coshtf = -r cosh(Z(c, T)). 

Comparing with (|5.3|) finishes the proof of the theorem. □ 
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There is no counterpart of Proposition Q] for spacelike vectors. Nevertheless, we obtain the following 
analogue of Lemma [6] . 



Lemma 7 Let c(s) be a spacelike solution of the Euler- Lagrange system (|5.1[) with A(s) 7^ 0. Then. 

1. The length of the velocity vector c(s) is constant along the solution; 

2. The horizontal coordinates are expressed by (|5.3|) . 

As the next step, we shall study the function A(s). First, let us prove some useful facts. 

Proposition 2 Let c(s) — (xi(s), X2(s), x 3 (s), Xi(s)) be a horizontal timelike (spacelike) curve. Then, 

1. -x\{s) - x 2 2 {s) + ±l(s) + ±l(s) = -a 2 (s)+ p 2 (s); 

2. c= a(s)T + b(s)X + uj(s)Y + w(s)N, with a = a, b = $, to = 0, w = a 2 - 1 . 

PROOF. Let us write the coordinates of c(s) in the basis T, X, Y, N as 

c(s) = a(s)T + (3(s)X + j(s)Y + S(s)N, 

where 





a = 


(c,T) 


= X 2 ±l — 


XxX-2 + X 4 X 3 - 


X 3 X 4 , 




P = 


(c,X) 


= 


— X4X2 + Xi± 3 


+ X2X4 


= 


1 = 


(c,Y) 


= X4X1 — 


X3X2 + X2X3 - 


XXX4, 


= 


5 = 


(c,N) 


= — XiX\ 


- X 2 X 2 + X 3 X 3 


+ X4X4 



By the direct calculation we get 

-a 2 + P 2 = -a 2 - S 2 + 1 + 7 2 = -x\ - x\ + x\ + x\. 
In order to prove the second statement of the proposition we calculate 

a — X2X1 — x\X2 + X4X3 — X3X4 = (c, T) = a, 

$ = —X3X1 — X4X2 + X1X3 + X2X4 = (c, X) = b. 
Differentiating the horizontality condition (|3.ip . we find 

= -^-(c, Y) = ^-{x4± 1 - X3X2 + x 2 x 3 - X1X4) = X4X1 - X3X2 + x 2 x 3 - X1X4 = (c, Y) = lu. 
Then, 

d d , \ 

= — (c, N) = -r{ — x\i,\ — X2X2 + x 3 x 3 + X4X4) = —x\X\ — X2X2 + x 3 x 3 + X4X4 

+{-x\ -xl+x 2 3 + x\) = (c, N) + {-a 2 + /3 2 ) = w - a 2 + /3 2 , 
by the first statement. The proof is finished. □ 

Theorem 4 The Lagrange multiplier A(s) is constant along the horizontal timelike (spacelike, lightlike) 
solution of the Euler- Lagrange system (|5.ip . 

PROOF. We consider the equivalent Lagrangian function L(x, x), changing the length function —a 2 + 
P 2 to — x\ — x\ + x\ + x\. The solutions of the Euler-Lagrange system for both Lagrangians give the 
same curve. Thus, the new Lagrangian is 

L(x, x) = - ( — x\ — x\ + ±1 + x 2 4 ) + A(s) (±1X4 — X4X1 — ±2X3 + x 3 x 2 ) ■ 
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The corresponding Euler-Lagrange system is 

— 'i\ — — AX4 — 2Ai4, 

— X2 = Ax 3 + 2Ai 3 , 
x a = —Xx2 - 2A± 2 , 

= — Xxi + 2A±i. 

We multiply the first equation by —x 4 , the second equation by X3, the third one by X2, and the last one 
by — x%, finally, sum them up. This yields 

£1X4 — £2^3 + £3X3 — £4X1 = X(x i + x 3 — x 2 — Xi) + 2\[XiXi + X3X3 — X2X2 — iixi) 

(c,Y) = -\ + 2\(c,N) => A = 0. 
We conclude that A is constant along the solution. □ 

We see from the proof of Lemma [6] that the function A(s) is just a linear function. This leads to the 
following property of horizontal timelike future-directed solutions of the Euler-Lagrange system (|5.1[) ■ 

Corollary 1 If c(s) is a horizontal timelike future- directed solution of (|5.ip . then the hyperbolic angle 
between its velocity and the time vector field T increases linearly in s. 

6. Hamiltonian formalism 

The sub-Laplacian, which is the sum of the squares of the horizontal vector fields plays the funda- 
mental role in sub-Riemannian geometry. The counterpart of the sub-Laplacian in the Lorentz setting 
is the operator 

C = i(-T 2 + X 2 ) = i ( - ( - x 2 d Xt + xA 2 + x 4 d X3 - x 3 d Xi f 

+ {x 3 d Xl + x 4 d X2 + xid X3 + x 2 d X4 ) 2 ^ . (6.1) 

Then the Hamiltonian function corresponding to the operator (|6.ip is 

1 / 2 2 

H{x, £) = - y - ( - x 2 £,i + xi& + £46 - £ 3 £ 4 ) + {X3£,l + + %l£,3 + X2U) 

= \{-r 2 + <?), (6.2) 

where we use the notations = d Xk , r ~ — a^Ci +xi^2 +^4^3 — a^d, and <; = X3^\ +0:4^2 +2:1^3 + 3:2^4. 
There are close relations between the solutions of the Euler-Lagrange equation and the solutions of the 
Hamiltonian system 

x - dH i - dH 

d£ ' x dx ' 

The solutions of the Euler-Lagrange system (|5.1|) coincide with the projection of the solutions of the 
Hamiltonian system onto the Riemannian manifold. In the sub-Riemannian case the solutions coincide, 
if and only if, the solution of the Euler-Lagrange system is a horizontal curve. We are interested in 
relations of the solutions of these two systems in our situation. The Hamiltonian system admits the 
form 

• dH 

x = — — = —TXJ + sxEi, 



i:)t 



n (6-3) 
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Lemma 8 The solution of the Hamiltonian system (|6.3p is a horizontal curve and 

t = a, q — (3, 
where a and (3 are given by (|3. 2|) and (|3 - 3[) respectively. 



(6.4) 



PROOF. Let c(s) = (xi(s), ^(s), %a(s), %4,[s)) be a solution of (|6.3p . In order to prove its horizontalitv 
we need to show that the inner product {x^xE^) vanishes. We substitute x from (|6.3|) and get 

(x, xE 2 ) = -t{xJ, xE 2 ) + q(xEi : xE 2 ) = 

by (gHH - 

Using the first line in the Hamiltonian system and the definitions of horizontal coordinates 
and (|3.3|) . we get 

a = (x, xJ) = —t(xJ 1 xJ) + q(xEx,xJ) = r, 
= (x, xE\) — —t(xJ, xEx) + s(xEi, xE\) = q 

from (PTU|) and (pOTj) . □ 



6.1. Geodesies with constant horizontal coordinates 

Lemma [8] implies the following form of the Hamiltonian system (I6.3[> 

xi = -a(-x 2 ) + (3x 3 , 
x 2 = -axi + f3x4, 

(6.5) 

±3 = — ax,4 + f3xi, 
x 4 = —a{—xz) + (ix2 1 

with constant a and (3. 
6.1.1. Timelike case 

In this section we are aimed at finding geodesies corresponding to the extremals (Section [5]) with 
constant horizontal coordinates a and (3 giving the vanishing value to the Lagrangian multiplier A. We 
give an explicit picture for the base point (1, 0, 0, 0). Left shifts transport it to any other point of AdS. 
Without lost of generality, let us assume that — a 2 + (3 2 = —1, a = cosh-0, (3 = sinh'0, where ip is a 
constant. 

The Hamiltonian system (|6.5[) written for constant a and j3 is reduced to a second-order differential 
equation 

ik = -Xk, k = 1, .. .4. (6.6) 

The general solution is given in the trigonometric basis as Xk = Ak cos s + Bk sin s. The initial condition 
a;(0) = (1,0,0,0) defines the coefficients Ak by A\ = 1, A 2 = A3 = A4 = 0. Returning back to the 
first-order system (|6.5p we calculate the coefficients Bk as B\ = 0, B 2 = —a, B 3 — (3, B 4 — 0. Finally, 
the solution is 

.Ti=coss, X2 =— cosh^ sins, X3 = sinh-i/'sins, 2:4 = 0. (6-7) 

These timelike geodesies are closed. Varying ip they sweep out the one-sheet hyperboloid xl+x 2 , — ! 2 = 1 
in M 3 . 

Let us calculate the vertical line T, the line corresponding to the vanishing horizontal velocity (a, (3) 
and with the constant value 7 = 1, passing the base point (1,0,0,0). Its parametric representation 
T = T(s) satisfies the system 
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a = X2X1 — X1X2 + X4X3 — = 0, 

(3 = —X3X1 — X4X2 + X1X3 + X2X4 = 0, 

7 = X4&I — X 3 X2 + X2X3 — X1X4 = 1, 

S — X\X\ + X2X2 — £3X3 — X4X4 = 0. 

The discriminant of this system calculated with respect to the derivatives as variables is (-1), and we 
reduce the system to a simple one 

Xi = -X4, x 2 =x 3 , x 3 =x 2 , x 4 = -xi, 

with the initial condition T(0) = x(0) = (1,0,0,0). The solution is 

r(s) = (coshs, 0, 0, — sinhs). 

The vertical line (hyperbola) T meets the surface (|6.7p at the point (1,0,0,0) orthogonally with respect 
to the scalar product in R 2,2 . Comparing this picture with the classical sub-Riemannian case of the 
Heisenberg group, we observe that in the Heisenberg case all straight line geodesies lie on the horizontal 
plane M 2 and the center is the third vertical axis. In our case the surface (|6.7p corresponds to the hori- 
zontal plane, timelike geodesies correspond to the straight line Heisenberg geodesies, and T corresponds 
to the vertical center. 

6.1.2. Spacelike/lightlike case 

Again we consider constant horizontal coordinates a and (3, and let us assume that — a 2 + 1 = 1, 
a = sinh-0, — cosh?/;, where tp is a constant. 

The Hamiltonian system (16.51) is reduced to the second-order differential equation 

x k =x k , fc = l,...4. (6.8) 

Arguing as in the previous case we deduce the solution passing the point (1,0,0,0) as 

xi = coshs, x 2 =— sinh-0 sinhs, x 3 = cosh -0 sinhs, X4 = 0. (6-9) 

These non-closed spacelike geodesies sweep the same hyperboloid of one sheet in K 3 . The vertical line 
r meets orthogonally each spacelike geodesic on this hyperboloid at the point (1,0,0,0). 

In the lightlikc case a 2 = 1 — 1 the Hamiltonian system (16. 5|) has a linear solution given by 

x\ = 1, X2 — —as, X3 — /3s, X4 = 0, 

which are two straight lines on the hyperboloid, and again T meets them orthogonally at the unique 
point (1,0,0,0). 

6.2. Geodesies with non-constant horizontal coordinates. 

If the horizontal coordinates are not constant, then we must solve the Hamiltonian system generated 
by the Hamiltonian (|6.2p . 

Fix the initial point x^ — (1, 0, 0, 0). We shall give two approaches to solve this Hamiltonian system 
based on a solution in Cartesian coordinates and on a parametrization of AdS. 

Solution in the Cartesian coordinates. It is convenient to introduce auxiliary phase functions 

Ul — Xl + X 2: U 2 =Xl-X 2 , U 3 =X 3 +X4, U4=X 3 -X4, 

and momenta 

Then the Hamiltonian (|6.2|) admits the form H = (— U41P2 + uiip 3 )(u 3 ipi — U21P4), and yields the Hamil- 
tonian system 
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Ul = U 3 (-U4lp2 +Uli) 3 ), Ul(0) = 1, 

U 2 = -U±(u 3 -ipi - U2lp4), U 2 (0) = 1, 

il3 = Ul(v,3lpl - m^i), U3(0) = 0, 

Ui = -U2{-Uii>2 + Ui0 3 ), w 4 (0) = 0, 



(6.10) 



for positions and 



(6.11) 



i>i = -^3(^3^1 - u 2 ip4:), V'i(o) = -4, 

02 = 1p4(-U4lp2 + "1^3)1 ^2(0) = B, 

■03 = -1pl(-U4^ 3 + Uil/l 3 ), 3 (O) = C, 
1p4 = Ip2(u34>l ~ U 2 1pi), 4>4(0) = D, 

for momenta with some real constants A,B,C, and D. For r and <; constant we get simple solutions 
mentioned in the previous section. We see that the system (|6.10f[6TTT|) has the first integrals 

"1^1 + «303 = A 
U 2 02 + «404 = B, 
"203 - "401 = C, 
"104 - W302 = D, 

and in addition, we normalize 0(0) so that the trajectories belong to AdS: u\u 2 + U3U4 = 1, and the 
Hamiltonian H = — 1 in the timelike case, in particular, the latter implies CD = 1. Then wc can deduce 
the momenta as 



01 


= Au 2 


- Cu 3 . 


"02 


= But 


— Du/t 


i>3 


= Cm 


+ Aui. 


■04 


= Du 2 


+ Bu 3 



Let us set the functions p = Ui/ux and q — u 3 /u 2 . Then substituting function in (|6.10|) . we get 

p = -{Dp 2 + (A- B)p+1/D), f»(0)=0, 

q = -{Cq 2 -{A- B)q + 1/C), q(0) = 0. 

The cases of the discriminant give the following options. Solving these equations for \A — B\ > 2, we 
obtain 

2 1 _ e - ay /(B-AF-i 

p{s) 



D 



(B-A- yftB - A) 2 - 4) — (B — A + \J (B - A) 2 - 4 ) e -V( s -^) 2 - 4 ' 



2D(1 - e-sy/tA-B)'-*) 

q{s) - 



(A-B- y/(A-B) 2 -A) -(A-B+ y/(A - B) 2 - 4)e _ V(^- B ) 2 - 4 
Next we use the relation u\ = —^114. Then, iii{pq + 1) = —pqux, and finally, 

Jo p(m{t) + 1 
f -MUZ 

p(t)q(t) + l 



U4 ( s )=p( s )exp / —EtMQ-dt. 



Taking into account u 2 = —ii 3 p, we get 



u 2 [s) = exp / eft, 
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, — q(t)p(t) 
u 3 (s)=a(s)e Xp ^ pmt) + 1 dt 



For A - B = 2 we get 



s -s 



Ui — (1 + s)e s , = (1 — s)e s , u 3 — —Dse s , 1*4 = — — e 
or in the original coordinates 

x\ = cosh s — s sinh s, £2 = — sinh s + s cosh s, 



2 V D J ' 2 V £> 

For A — B = —2 and for |A— i?| < 2 one obtains the solution analogously in the timelike case CI? = 1. 
Thus we get a two-parameter D and A — B family of geodesies passing through the point (1,0,0,0). 
The parameters D and A — B have a clear dynamical meaning. Namely, 

D = -u 3 (0) = -(± 3 (0) + ±4(0)), C = i = -«4(0) = -(2:3(0) - ± 4 (0)), 

and 

u 3 (0) u*(0) i 3 (0) + i 4 (0) x 3 (0)-x A (0)' 
The spacelike case CD = —1 is treated in a similar way, but we omit awkward formulas. 

Parametric solution. We present the parametric form of timelike and spacelike geodesies starting 
from the point (1, 0, 0, 0). The forms of solutions with constant velocity coordinates (|6.7p and (|6.9p give 
us an idea of a suitable paramctrization for geodesies with different causality. 
Timelike geodesies. 

We use the parametrization in a neighborhood of (1, 0, 0, 0), given by 

xi = cos 4> cosh xi, 
x 2 = sin cosh X2, 

(0.12) 

x 3 = sin </> sinh X2, 
x 4 — cos sinh , 

where 4> £ (— ? , J), XI1X2 € (00,00). We note that the timelike solution with constant velocity coor- 
dinates (|6.7p followed from this parametrization if we set <fi = — s, xi = 0j an d X2 = ~~ "0- The vertical 
line r is obtained by setting <j) = 0, Xl = — s > an d X2 = 0. 

In this parametrization the vector fields T, X, and K admit the form 

T = 2cosh(xi - X2)d<i > + d Xl tan0sinh(xi - X2) + d X2 cotan</>sinh(xi - X2), 
X = 2sinh(xi - Xa)^ + d Xl tan0cosh(xi - X2) + d X2 cotan</>cosh(xi - X2), 
Y = d Xl -d X2 . 

The vertical direction is given by the constant vector field Y. Let c(s) — (4>(s), x( s ), X2( s )) be a 
curve starting at c(0) = (0, 0, X2(0)). The horizontal coordinates (|3.2|) and f|3.3|) with respect to given 
parametrization are 

a = -0cosh(xi - X2) + ^(Xi + X2) sin(20) sinh(xi - X2), 
(3 = 0sinh(xi - X2) + 7^X1 + X2) sin(20) cosh(xi - Xi)- 
Then, the square of the velocity vector c(s) is 



~a 2 +P 2 = -4? + \(xi + X2) 2 sin 2 (20). 
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The speed is preserved along the geodesies and is equal to the initial value at the point (1, 0, 0, 0), or in 
our parametrization (0,0,X2(0). Therefore, 

(c(O),c(O)) = (- a 2 +/3 2 )(O) = -0^(O), 

and we obtain timelike geodesies starting from (0,0,X2(0)) if 0(0) ^ 0, and lightlike geodesies in the 
limiting case 0(0) = 0. 

The Hamiltonian H associated with the operator 



becomes 



C = \{-T 2 + X 2 ) = l(-4d 2 + tan 2 0£^ + cotan 2 0d 2 2 + 2d Xl d X2 ), 
H(<t>, Xu X2, V, &» 6) = \(~^ 2 + & tan 2 + £f cotan 2 + 2£i£ 2 ) , 



where we set = ip, d Xl = £i, and d X2 
The Hamiltonian system 



Xi = £,i tan 2 <j) + £, 2 , 
X2 = 6 cotan 2 + 6, 
= -4V> 3 

6 = o, 

6 = 0: 



(6.13) 



. 9 tan0 2 cotan 

^ = 2~I + £2 • 2 7 ■ 

cos 2 snr 

shows that £1 and £2 are constants. If both constants vanish, then we get 

*i = 0, X2 = 0, 0=-4V>, = 0, 

which leads to the trivial solution (|6.7[) . Since we are looking for a solution in a neighborhood of 
(0, 0,%2(0)), we put £2 = 0. Let us solve the Hamiltonian system (|6.13[) with the initial conditions 

0(0) = 0, xi(o) = o, X2(0) = x 2 °\ ^(o)-0 (o) , 6(0) = a, 6(0) = 0. 

From the third and from the last equations we get = —4-0 = 4£ 2 ^r^. Multiplying by and integrating 
we obtain 

2 (s) = C 2 + 4£ 2 tan 2 0(s), C = 2 (O) = 16?/> 2 (0). (6.14) 
Let us assume C 2 > 0. Simplifying (|6.14[) . we come to 

cos d<t> 

= ±ds. 



v /c 2 + (4^ 2 -C* 2 )sin 2 



According to the sign of 4£ 2 — C 2 , one gets three different types of solutions. 

Case 1: 4£ 2 — C 2 = 0. Integrating from to some value of s we get the solution in the form sine 

±|C|«. 

Case 2: 4£ 2 - C 2 > 0. The solution follows as 



4£ 2 - C 2 



C 2 

Case 3: 4£ 2 — C 2 < 0. The solution is obtained as 



sin0 = ±sinh(s v /4e 2 - C 2 ). 



° 2 c ^ sin0 = ±sin( SV / C 2 - 4£ 2 ). 

In order to calculate the value of Xii we express tan 2 from the Cases 1-3 and integrate the first 
equation of the Hamiltonian system. Observe that xi — £1 is constant and 0(0) = — 4-0' -' 7^ 0. The 
following theorem is proved. 
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Theorem 5 The timelike geodesies starting from the point 4>(0) = 0, Xi(0) = 0, X2(0) = xi w ith 
some <f>(0), a constant value of \2, and an arbitrary Xi( s ) satisfy the following equations: 

if^xl = 2 (O) then 



sin0(s) = ±|C|s, 

' i+0(g)s 

l-0(O)s 

(0) 



Xi(s) = -X2S + ^yln 



• Xa(s) =X2S + x 2 / 
z/4xi > </> 2 (0) tfien 

• sin d>(s) = ± 



• Yi fsl = -Y9S + Y9 f" 4 xI-^(°) ds 

^ • (o) 

• X2(s) = x2s + x 2 ; 

and if f 4x1 < 4> 2 {0) then 

• sin </>fsl = ± 



^/~^ si „( sv ^T^), 



•^"-»-^ ^,.,„;;;;:,». l -,» i ^ 

C ^ • i (0) 
• X2(S) =X2S + X2 ■ 

The integrals can be easily calculated and they involve trigonometric and hyperbolic functions, and 
depend on the relations between 4x|, <^> 2 (0). 

Spacelike geodesies. 

We use another parametrizaition in a neighborhood of (1, 0, 0, 0) suitable in this case 

x\ = cosh cosh xi, 

x 2 = sinh0coshx2, ,„ ir . 

(6.15) 

£3 = sinh0sinhx2, 
Xi = cosh <f> sinh xi , 

where 4>, xi, X2 € ( — °°j 00). Observe that the spacelike solution with constant velocity coordinates (|6.9[) 
follows from this parametrization if we set 4> = s, xi = and X2 = — "0- The vertical line T is obtained 
as previously, by setting <f> = 0, Xi = ~ s ; an d X2 = 0. 
The vector fields T, X, and Y become 



T = 2sinh(xi - X2)d < f> - d Xl tan0cosh(xi - X2) + d X2 cotan <j> cosh(xi - X2), 
X = 2cosh(xi - X2)d<t, - d Xl tan0sinh(xi - X2) + d X2 cotan </>sinh(xi - X2), 
Y = d Xl -d X2 . 

The vertical direction is again given by a constant vector field Y, Let c(s) = (0(s), x(s), X2(s)) be a 
curve such that c(0) = (0,0, X2(0)). The horizontal coordinates (|3.2|) and (|3.3j) with respect to this 
parametrizaition are 

a = 0sinh(xi - X2) - \(xi + X2) sinh(20) cosh(xi - X2), 

i 

(3 = 0cosh(xi - X2) - g (Xi + X2) smh(20) smh(xi - X2)- 
Then the square of the velocity vector c is 

-a 2 + /3 2 = 2 - + X2) 2 sinh 2 (20). 
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Since the speed is preserved along the geodesies, it is equal to 2 (O), and we obtain spacelike geodesies 
starting from (0,0,X2(0)) for 0(0) + 0. 

The Hamiltonian H associated with the operator 

1 ( _ T 2 + X 2 )= 1 ( 

becomes 



C = -(-T 2 + X 2 ) = -(Ad 2 - tanh 2 4>d 2 Xl - cotanh 2 0d 2 2 + 2d Xl d X2 ) 



H{cf), xi, X2>, Ci. 6) = ^(4^ 2 - e 2 tan 2 - g cotan 2 + 2£i&), 



where we set 3^ = V*> ^xi = £l> an d ^X2 = £2- 
As in the previous case, the Hamiltonian system 





= — £i tanh 2 


X.2 


= — £2 cotanh 2 




= 


ii 


= o, 




= o, 




^ 2 tanh 


1 cosh 2 



(6.16) 



^ 2 cotanh 
2 sinh 2 

gives that £i and £2 are constants. If both constants vanish, we get 

Xi =0, X 2 = 0, = -4^>, </> = 

which leads to the spacelike trivial solution. Setting £2 = 0, we solve the Hamiltonian system (|6.16[) 
with the initial conditions 

0(0) =0, xi(0) = 0, X2 ( )= x ( 0) , V^(o) = v (0) , £i(0) = fi, 6(0) =0. 

An analogue of (|6.14p is 

2 (s) = C 2 + 4£ 2 tanh 2 0(s), C = 2 (O) = 16^ 2 (0) ^ 0. (6.17) 

Arguing as in the timelike case, we prove the following statement. 

Theorem 6 The spacelike geodesies starting from the point 0(0) = 0, Xi(0) = 0. X2(0) = X2 ™th 
some 0(0), a constant value of X2> and an arbitrary Xi( s ) have the following equations: 



Sillh<KS) = ± V 2 (O) ( + O) 4xl Sin ^ S V^(°) + 4 ^ 



cxi(s) = ~X2S + arccotanh ^ 4 * 2 cotan (sy / 2 (O) + 4x1)) > 

X2(s) =X2S + X 2 0) - 



7. Geodesies with respect to the distribution T> — span{X, Y} 



This case reveals the sub-Riemannian nature of such a distribution. In principle, one can easily modify 
the classical results from sub-Riemannian geometry (Chow-Rashevskii theorem, in particular). However 
we prefer to modify our own results proved in previous sections to show some particular features and 
to compare with the sub-Lorentzian case defined by the distribution T> = span{T, X}. 

Lemma 9 A curve c(s) = (xi (s), X2(s), Xa(s), x±(s)) is horizontal with respect to the distribution T> = 
span{A, Y}, if and only if 

X2X1 — X1X2 + — X3X4 = or (xJ,c)=0. (7-1) 
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PROOF. The tangent vector to a curve c(s) = (xi(s),X2(s),x^(s),Xi(s)) written in the left-invariant 
basis is of the form 

c(s) = aT + /3X + iY. 

Then 

a = (c,T) — Tc ■ T — x\X2 — X2X\ + X3X4 — ±4X3. 
We conclude that a = 0, if and only if, (|7.1[) holds. □ 

In this case a curve is horizontal, if and only if, its velocity vector is orthogonal to the vector field T. 
The left-invariant coordinates f3(s) and "/(s) of a horizontal curve c(s) = (xi(s), X2(s), Xa(s),Xi(s)) are 

(3 = (c, X) = -2:3X1 — + X1X3 + x 2 X4 = (xEi,c). (7.2) 

7 = (c, y) = -X4X1 + x 3 x 2 - x 2 x 3 + X1X4 = (xE 2 ,c). (7.3) 

The form w = —X2<ix\ + X\dx2 — x^dx^ + x 3 dx^ — ~(xJ,dx) is a contact form for the horizontal 
distribution V = span{X, Y}. Indeed 

w(N) = 0, w(T) = 1, w(X) = 0, w(Y) = 0. 

Thus, kerw = spanjiV, X, Y}, The horizontal distribution can be defined as follows 

V = {V e TAdS : w(V) = 0}, or V = keru; n TAdS. 

The length 1(c) of a horizontal curve c(s) : [0, 1] — > AdS is given by 

1(c) = [ {c(s),c(s)) 1 / 2 ds= [ {0 2 (s) +1 2 ( s )) 1/2 ds. 
Jo Jo 

The restriction of the non-degenerate metric (•, •) onto the horizontal distribution T> C TAdS gives a 
positive-definite metric that we still denote by (•, -)x>. Thus from now on, we shall work only with one 
type of the curves (that we shall call simply horizontal curves) , since the horizontality condition requires 
the vanishing coordinate function of the vector field T. 

7.1. Existence of horizontal curves 

The following theorem is an analogue to Theorem Q] proved for the distribution T> — span{T, X} in 
Section 4. 

Theorem 7 Let P, Q € AdS be arbitrary given points. Then there is a smooth horizontal curve 
connecting P with Q. 

PROOF. We use parametrisation (|4.ip . in which the horizontality condition for a curve c(s) is ex- 
pressed by (|4.3p as 

ip + <p cosh 20 = 0. 
This equation is to be sold for the initial conditions 

c(0) = P or <p(0)=tp Q , #))=^o, 0(0) = 6 a , 

c(l) = Q, or (p(l)=<pi, V(l) = ^i, 9(1) =0i. 

Let ip = tp(s) be a smooth arbitrary function with ijj(0) — lim i/^s) and "0(1) = li m ^(s)- Set 

26*(s) = arccoshp(s). Then the equation ()4.3j) admits the form 

- tp(s)=- — __ + (/ 5(0). 



cosh 26 p(s) ' 7 p(s) 
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Denote q(s) = Since q(0) = ^^"^ , = ■ and Jo ^( s ) ds = <Po - ¥>i applying Lemma E] 

we conclude that there exists such a smooth function q(s). The function p(s) is found as p(s) = -^jy- 
We get a curve c(s) = ((p(s),ip(s),9(s)) with 

ip = i/j(s), 



i 

9(s) = — arccoshp(s) 



□ 



Remark 3 Observe that in the general Chow-Rashevskii theorem smoothness was not concluded. 

ih 

<Po—<Pi ' 



Theorem 8 Given two arbitrary points P = P((fo, ipo, 9q) and Q — $o) with 29q = arccosh — ^ 



there is a horizontal curve with the constant 9-coordinate connecting P with Q. 

PROOF. If the ^-coordinate is constant, then the governing equation is 

ip = —if cosh 29q ip(s) — — f{s) cosh 26*o + C. 

Applying the initial conditions 

c(0) = (ip Q ,ipo,6o), and c(l) = (ipi, ipi, 6 Q ) , 

we find 

2#o = arccosh I , G = Vo + Vo • 

Therefore, for any parameter ip, the horizontal curve 



29n — arccosh . 



(f - ip i / tf - (f X 

joins the points P = P(ifo,ipo,9 ) and Q — Q((pi,ipi,0o). □ 

7.2. Lagrangian formalism 

Dealing with T> = span{X, Y} and a positive-definite metric (•, -)x> on it, one might compare with 
the geometry generated by the sub-Riemannian distribution on sphere S 3 in [3] . The minimizing length 
curve can be found by minimizing the action integral 



2 jo 

under the non-holonomic constrain a — (c, xJ) = 0. The corresponding Lagrangian is 

L(c,c) = \{P 2 {s) + 7 2 (s)) + X(s)a(s). (7.4) 

The extremal curve is given by the solution of the Euler-Lagrange system (|5.ip with the Lagrangian (|7.4|l . 

Let us make some preparatory calculations. Write the system (15. ip for the Lagrangian (|7.4p as the 
follows 

2/3±3 + 27x4 — 2A±2 + /3x3 + 7x4 - Ax2 = 0, 

2/3x4 — 27x3 + 2A±i + /3x4 — 7x3 + Axi = 0, 
— 2/3xi + 27x2 — 2AX4 — 0x\ + 7x2 — AX4 = 0, 
—2/3x2 — 27x1 + 2A±3 — /3x2 — 7x1 + AX3 = 0. 
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Multiply the equations by X3, Xi,xi, and X2, respectively and sum them up. We get 

2/3{c,N) - 2 7 (c,T) - 2\{c,Y) - /3 + O7 + OA = => (i = 2A 7 , 
2/3(c,T) - 2 7 (c, AO +2A(c,X) + 0/3- 7 + 0A = => 7 = 2A/3. 

Let us consider two cases. 

Case A(s) = 0. In this case equation (|7.2I) admits the form 

$ = 0, 7 = 0, (7.5) 

and we deduce the following theorem. 

Theorem 9 There are horizontal geodesies with the following properties: 

1. The coordinates a = (c, T) = 0, /3 = (c, X), and 7 = (c, Y) are constant; 

2. The length |c| along the geodesies; 

3. The angles between the velocity vector and horizontal frame is constant along along the geodesic. 



PROOF. Taking into account the solution of (|7.5[) . we denote f3(s) = (3 and 7 (s) = 7. Then the length 
of the velocity vector |c| = \J ft 2 + 7 2 is constant. 

Since (c, X) = (c,X)x> — \c\d\X\d cos(Zc, X), (c,Y) = (c,Y)x> = \c\v\Y\v cos(Zc, F), we have 

cos(Zc,X)= cos(Zc,T) 



that proves the third assertion. □ 



Case A(s) 7^ 0. 

Theorem 10 There are horizontal geodesies with the following properties: 

1. The velocity vector \c\ of a geodesic is constant along the geodesic; 

2. The angles between the velocity vector and the horizontal frame are given by 

Zc,X = cs + 9 , Zc,Y = ^ - cs + O . 



PROOF. Since 

f3 = 2A7, 7 = 2A/3 (7.6) 

implies Jj(/3 2 + 7 2 ) = 0, we conclude, that the length of the velocity vector \c\ is constant. Taking into 
account positivity of (3 2 + 7 2 let us denote it by r 2 . Set (3 = rcosO(s) and 7 = rsin#(s). Substituting 
them in (|7.6p . we get 

9(s) = 2A(s) => 9{s) = 2 J A(s) ds + 9 - 

Let us find the function A(s). Observe that 

f3 2 + 7 2 = -i 2 - x\ + i 2 + i 2 . 

It can be shown similarly to the proof of Proposition [21 having a = 5 = 0. By the direct calculation 
(see also Proposition ^ we show that 

(c,T) = ±(c,T) = 0. 
Now, we consider an equivalent to (17. 4| extremal problem with the Lagrangian 

L(c,c) = i(-i 2 -i 2 + i 2 +i 2 ) +A( S )(c,T). (7.7) 
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The Euler-Lagrange system admits the form 

— X\ = —2\X2 — A.T2, 

— X2 = 2Xii + Xxi, 

X3 = — 2A±4 — \X4, 
X4 = 2A±3 + A1E3. 

Multiplying these equations by x%, —x±, —X4, X3 respectively and then, summing them up, we obtain 

-<c,T) = 2X(c,N) - A. 

This allows us to conclude, that the function A(s) is constant along the solution of the Euler-Lagrange 
equation that yields the second assertion of the theorem. □ 

7.3. Hamiltonian formalism 

The sub-Laplacian is £ = X 2 + Y 2 and the corresponding Hamiltonian function is 

1 / 2 2\ i 

H(Z,Q = X \i x 3^1 + X4& + X%^3 + X2U) + ~ %3& ~ X 2 £ 3 + Xt&) j = - {<? + K 2 ). 

The Hamiltonian system is written as 



(7. 



dH 

x = —r = qxEi + kxE 2 , 
ox 

As in the previous section we are able to prove the following proposition. 
Proposition 3 The solution of the Hamiltonian system is a horizontal curve and 

? = (3, k = 7. 

Corollary 2 The Hamiltonian function is the energy H(x,^) = \{(3 2 + 7 2 )- 
7.4. Geodesies with constant horizontal coordinates. 

In this section we consider constant horizontal coordinates (3 and 7. Making use of Proposition [3] we 
write the first line of the Hamiltonian system (|7.8|) in the form. 

xi = (3x 3 + 7x4, 

X 2 = (3X 4 ~ 7X3, 

(7.9) 

X 3 = (3X4 - 7X2, 



Xi = (3x 2 + ix 



We give an explicit picture for the base point (1, 0, 0, 0). Without lost of generality, let us assume that 
/3 2 +7 2 = 1, (3 = cos?/;, 7 = simp, where ip is a constant. 

The Hamiltonian system (|7.9p written for constant (3 and 7 is reduced to a second-order differential 
equation 

x k =x k , k=l,..A. (7.10) 

The general solution is given in the hyperbolic basis as Xk = A k cosh s + B k sinh s. The initial condition 
x(0) = (1,0,0,0) defines the coefficients A% by A\ = 1, A% = A3 = Ah = 0. Returning back to the 
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first-order system (|7.9I) we calculate the coefficients Bk as B\ = 0, B2 = 0, B3 = (3, B4 = 7. Finally, 
the solution is 

xi = coshs, X2 = 0, X3 = cos"0sinhs, 2:4 = sin^sinhs. (7-11) 

Varying ip they sweep out the two-sheet hyperboloid X-y X^ X ^ = 1 in R 3 . We use only one sheet 
containing the point (1,0,0,0). Geodesies are hyperbolas passing this point. 

The vertical line corresponds to the vanishing horizontal velocity (/3, 7) and with the constant value 
a = 1, passing the base point (1, 0, 0, 0). The solution is 

T(s) = (cos s, sins, 0,0). 

The vertical line (circle) T meets the surface (|7.11j) at the point (1,0,0,0) orthogonally with respect to 
the scalar product in R 2 ' 2 . 

7.5. Geodesies with non- constant horizontal coordinates. 

If the horizontal coordinates are not constant, then we must solve the Hamiltonian system generated 
by the above Hamiltonian. 

Solution in the Cartesian coordinates. Fix the initial point x^ = (1,0,0,0). In the Cartesian 
case it is convenient to introduce complex coordinates z = x\ + 1x2, w = X3 + ix<±, ip = £1 + i£,2, and 
ip = £3 + ^4. Hence, the Hamiltonian admits the form H = \zip + wip\ 2 . The corresponding Hamiltonian 
system becomes 

z = w{ztj} + wep), z(0) = 1, 
w = z{zijj + w(p), w(0) = 0, 
tp = —ijj(zijj + w(p), (p(0) — A — iB, 
ip = —(p{ztp + wip), V'(O) = C — iD. 

Here the constants A, B, C, and D have the following dynamical meaning: w(0) = C + iD, and 2B = 
iw(0)/w(0). This complex Hamiltonian system has the first integrals 

zip + wip — C + iD, 
zip + wip — A — iB, 

and we have \z\ 2 — \w\ 2 = 1 and H = C 2 + D 2 = 1 as an additional normalization. Therefore, 

<p = z(A + iB) - w(C + iD), 

ip = z(C + iD) - w(A + iB). 

Let us introduce an auxiliary function p = w/z. Then substituting ip and ip in the Hamiltonian system 
we get 



p(s) = -(C-iD)- 



Vl-B 2 -iB + Wl-B 2 - iB)e- 2s ^^ 
Taking into account that zz = ww, we get the solution for B ^ 1 

f P(t)p(t) , 
Z(S) = 6XP Jo l-b(*)l 2 ' 

and 

For B = 1 the solution is 

z(s) = (1 + is)e ls , w{s) = s(C + iD)e~ ls . 
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Parametric solution. Let us present the parametric solution in this case. We use the parametrization 
in a neighbourhood of (1, 0, 0, 0) given by 

Xx = cosxi cosh0, 
x-> = sin Yi cosh0, 

(7.12) 

X3 = cos X2 sinh <f>, 
.T4 = sin xi sinh (f>, 

where <j> S (—00,00), XI7X2 S (— §, §)■ We observe that the solution with constant velocity coordi- 
nates (|7.11| follows from this parameterization when we set <j> = s, Xi = 0, and xi = 0- The vertical 
line (circle) T is obtained by setting = 0, xi = s , an d X2 = 0. 

In this parametrization, the vector fields T, X, and Y admit the form 

T = 8 - d 

X = 2cos(xi - X2)^ - d Xl tanh0sin(xi - X2) + d X2 cotanh </>sin(xi - X2), 

Y = 2sin(xi - X2)d cf) - d Xl tanh0cos(xi - X2) + d X2 cotanh </>cos(xi - X2). 

The vertical direction is given by the constant vector field T. 
The Hamiltonian H associated with the operator 

C=^(X 2 + Y 2 ) = ~(4d| + tanh 2 4>d\ x + cotanh 2 <\>d\ 2 - 2d Xl 8 X2 ) 

is given as 

H(4>, Xi, X2, V, 6, &) = ^(4^ 2 + C 2 tanh 2 + £ 2 2 cotanh 2 - 2&£ 2 ), 

where we set = ip, d Xl = £1, and d X2 = £2- 

Description of geodesies is collected in the following theorem. 

Theorem 11 The geodesies starting from the point (f)(0) = 0, Xi(0) = 0, X2(0) = X2 with some 0(0), 
a constant value of X2, and an arbitrary Xi(s) have the following equations. 

#4x! = 2 (O) then 

• sinh0(s) = ±|C|s, 

• Xi(s) = X2S - arctan 0(O)s, 

1 \ 1 (°) 

• X2(S) = -X2S + X2 ■ 

If 4x1 > ^(O) then 



. s inh0( S )=±J i -|^sin( Sv /4 X 2 -^ ( o) 



• Xi 00 = X2S - 5^- arctan (^mE|M tan (^x 2 _ ^))) , 

i \ (0) 

• X2(SJ = -X2S + X2 • 

#4x1 < 2 (O) tten 



sinh </,(*) = ±\\ ^% sinh ( n /^ 2 (0) - 4x 



t2 J; 



• Xi 00 = X2.s - ^ arctan ( ^"ff* 3 cotan (^4x1-^(0))), 

f ^ 1 (°) 

• X2(«) = -X2S + X2 • 
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